Scaling anomalies in the coarsening dynamics of fractal viscous fingering patterns 
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We analyze a recent experiment of Sliaron et al. (2003) on the coarsening, due to surface ten- 
sion, of fractal viscous fingering patterns (FVFPs) grown in a radial Hele-Shaw cell. We argue that 
an unforced Hele-Shaw model, a natural model for that experiment, belongs to the same univer- 
sality class as model B of phase ordering. Two series of numerical simulations with model B are 
performed, with the FVFPs grown in the experiment, and with Diffusion Limited Aggregates, as 
the initial conditions. We observed Lifshitz-Slyozov scaling t^^^ at intermediate distances and very 
slow convergence to this scaling at small distances. Dynamic scale invariance breaks down at large 
distances. 
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Coarsening is an important paradigm of emergence of 
order from disorder. It has been extensively studied in 
two-phase systems quenched from a disordered state into 
a region of phase coexistence Q, |^ |^ . In another class 
of systems, disordered configurations are generated by 
an instability of growth in combination with noise, and 
they often exhibit long-range correlations and fractal ge- 
ometry 01 . Examples include fractal clusters developing 
in the process of solidification from an under-cooled liq- 
uid (a, fractal clusters on a substrate grown by deposi- 
tion jg and fractal viscous fingering patterns (FVFPs) 
formed by the Saffman- Taylor instability in the radial 
Hele-Shaw cell 0- When the driving stops, the frac- 
tal clusters coarsen by surface tension, and the coarsen- 
ing dynamics provide a valuable characterization of these 
systems. 

An important simplifying factor in the analysis of 
coarsening dynamics is dynamic scale invariance (DSI): 
the presence of a single, time-dependent length scale L{t), 
so that a normalized pair correlation function C (r, t) de- 
pends, at long times, only on r/L(t). The coarsening 
length scale L[t) often exhibits a power law in time 
For systems with short-range correlations there is a lot 
of evidence, from experiments and numerical simulations, 
in favor of DSI 3]. For systems with long-range corre- 
lations the situation is more complicated. In the case 
of a non-conserved order parameter DSI was established 
in particle simulations following a quench from T — Tc 
to T = 0- Implications of mass conservation in DSI 
were addressed more recently, in the context of coarsen- 
ing of fractal clusters. Most remarkable of them is the 
predicted decrease of the cluster radius with time . As 
of present, only the systems where the conservation law is 
imposed globally, rather than locally, have been found to 
indeed show this effect . On the contrary, the "frozen" 
structure of fractal clusters at large distances, observed 
in simulations of locally conserved (diffusion-controlled) 
fractal coarsening [iH |l^ implies breakdown of 

DSI in these systems [ill The frozen structure is 



due to Laplacian screening of transport at large distances 
0. 

An additional scaling anomaly, observed in the numer- 
ical simulations of diffusion-controlled fractal coarsening 
Ull , was the presence of two apparently dif- 

ferent dynamic length scales. For one of them, deter- 
mined from the time-dependence of either the slope of 
the Porod-law part of C(r) [H E 13, or the cluster 
perimeter ll0llllU l3|. a power law in time was reported: 
Li ^ ^0. 20-0. 23 Another length scale, determined from a 
knee-like structure in C (r) at moderate distances behaves 
like La - J^. While L2(i) can be identified as 

Lifshitz-Slyozov length scale ~ i^/^ the length scale 
Li{t) looks unusual. 

Strikingly similar results were recently obtained in ex- 
periment on the coarsening dynamics of a different sys- 
tem: radially grown FVFPs in a Hele-Shaw cell The 
frozen structure at large distances, observed in Ref. [l^ . 
clearly indicates breakdown of DSI. Furthermore, two dif- 
ferent time-dependent length scales, with apparent dy- 
namic exponents 0.22 and 0.31, are observed Why 
is this system so similar to the diffusion-controlled sys- 
tem? Where does the exponent 0.20 — 0.23 come from? 
These questions are addressed in the present work. We 
first suggest an unforced Hele-Shaw model and discuss its 
properties. A scaling argument indicates that this model 
belongs to the same universality class as the so called 
model B, the standard model of the diffusion-controlled 
phase separation in two-phase systems ■ Assuming uni- 
versality, we performed a series of numerical simulations 
with model B, where the FVFPs, grown in experiment of 
Sharon et al. jl^, are used as the initial conditions for 
the minority phase. Then we report additional simula- 
tions of fractal coarsening, with DLAs (Diffusion Limited 
Aggregates) as the initial conditions. These two series of 
simulations show Lifshitz-Slyozov scaling t^^^ at interme- 
diate distances. Breakdown of dynamic scale invariance 
at large distances is confirmed. However, the existence 
of an anomalous power law in Li(t) is disproved. 
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A natural description of coarsening of FVFPs is pro- 
vided in terms of an unforced Hele-Shaw (UHS) flow. 
Consider a Hele-Shaw flow and assume that the driv- 
ing fluid (for example, air) has negligible viscosity, so 
that the pressure inside it is spatially uniform. When 
the plate spacing b is very small, the flow is effectively 
two-dimensional, and the velocity of the viscous fluid (for 
example, oil) is v (r, t) — — (6^/12/i) Vp (r, t), where p is 
the pressure and /i is the dynamic viscosity of the driven 
fluid. As the fluids are immiscible, the interface speed is 



I2n 



(1) 



where index n denotes the components of the vectors 
normal to the interface, and the gradient is evaluated 
at the respective points of the interface 7. Assuming 
incompressibility of the driven fluid, V- v = 0, one arrives 
at Laplace's equation for the pressure: 



V^p = . 

The pressure jump across the interface is |l6j 



cr 



1 + 3.8 



(^) 



2/3' 



(2) 



(3) 



where a is surface tension, and K is the curvature of the 
interface. At the coarsening stage the interface speed is 
very small, so the second term in the square brackets 
can be neglected. The first term does not depend on the 
coordinates, so one arrives at a Gibbs-Thomson relation 



Ap = - cr/C . 



(4) 



To close this set of equations, one more condition is 
needed. A natural condition to demand during the 
growth stage is a constant- in-time driving pressure [l7j . 
or a constant areal flow rate of the driving fluid. Each of 
these conditions assumes evacuation of the driven fluid 
at the external boundary of the system. In the coars- 
ening problem both the supply of the driving fluid, and 
evacuation of the driven fluid are blocked. Therefore, the 
normal component of the velocity of the driven fluid at 
the external boundary F should vanish, which follows 



V„p| 



0. 



(5) 



Equations O, @, (0)1 and © define a one-sided version 
of the UHS model. Similar models have been used in 
the context of break-ups (pinc h-offs) of bubbles, driven 
only by surface tension llSl ■ The UHS model has two 
important properties: (i) The total area A of the driving 
fluid is constant, (ii) The total length of the interface is 
a non- increasing function of time [20j . 

Now let us compare the UHS model with model B, the 
phase-field formulation of which is given by the Cahn- 
Hilliard equation for the order parameter w(r,t) Q: 

^^+\V'{V''u + u-u^)^Q. (6) 



At late times, the two-phase dynamics are describable by 
an asymptotic sharp- interface theory |2lj |. In the sharp- 
interface limit, the interface speed is 



(7) 



where potential $(r,i) is a harmonic function in each of 
the two phases in and out. The boundary conditions are 
<^\^ = (%/2/3) JC and V„$ |r = 0. 

How are these two problems related? To begin with, 
the sharp-interface limit of model B has the same prop- 
erties (i) and (ii) as the UHS model |21jj, so each of the 
two models describes interface-shortening dynamics un- 
der area conservation. The models do differ from each 
other considerably in the final outcomes of the coarsen- 
ing dynamics. For a steady state solution of the UHS 
model one has simply p = const. Therefore, possible 
stationary shapes of domains of the driving fluid in the 
UHS model are one or more circular bubbles of arbitrary 
radii. On the contrary, in model B, $ = const cannot be 
a steady state solution in the presence of more than one 
bubbles, because it cannot obey all the boundary condi- 
tions on the multiple interfaces. Therefore, a generic final 
state here is always a single circular bubble. In model B 
bubbles compete for material via diffusion through the 
majority phase. Obviously, this competition mechanism 
(Ostwald ripening 0,01) is absent in the UHS model. 

This difference between the two models becomes cru- 
cial after the driving fluid breaks up into multiple bub- 
bles. Before it happens, the two models can be expected 
to behave similarly. A simple argument for this follows 
from scaling analysis. Consider coarsening of a domain 
of complex shape and assume for a moment DSI, that 
is a single relevant length scale L = L(t). The pressure 
jump across the interface can be estimated from Eq. Q): 
Ap - a/L. Then, from Eq. (P), w„ ~ b^a/{iiL'^). On 
the other hand, u„ ~ L. This yields a coarsening law 
L{t) ~ {b^at/ fiY^^ . We checked that this estimate is in 
excellent agreement with the experimental result |l4| for 
L2{t), for two latest decades of time. 

The same power law t^^^ is obtained in model B 0,0]. 
Therefore, if DSI holds, the two models belong to the 
same universality class. In reality, each of these two sys- 
tems exhibits breakdown of DSI at large distances, when 
one deals with fractal clusters at ^ = [Tllllllli. . 
At intermediate distances, however, the classic exponent 
1/3 is observed in both systems |l3l Il4l|. Therefore, we 
conjecture that, prior to major breakup, the two models 
belong to the same universality class. 

Based on this conjecture, we performed two series of 
simulations with model B [Eq. ©J. Details of our 
numerical procedure and diagnostics can be found in 
Refs. hUhJ. In the first series of simulations we used 
FVFPs grown in experiment 0| as the initial conditions 
for the "minority phase" u = 1. The fractal dimension 
of these patterns, determined from the pair correlation 
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function, is close to 1.71. The (scaled) system size was 
1024 X 1024, with periodic boundary conditions. The 
(scaled) time range of the simulations was < t < 3- lO"'. 
Figure 1 shows snapshots of the simulated coarsening dy- 
namics. The snapshots closely resemble those observed 
in experiment 14J. Figure 2 presents the (normalized) 
equal-time pair correlation function C (r, t) at different 
times, and the characteristic dynamic length scales. The 
data is averaged over 7 simulations with different FVFPs. 
C(r, t) in a linear scale is shown in Fig. 2a. At small dis- 
tances C goes down linearly with r (the Porod law) |3|, 
and the inverse slope of this linear dependence yields the 
"coarsening length scale" Li{t) depicted in Fig. 2d. The 
log-log plots of C{r,t) (Fig. 2b) indicate an invariable 
fractal dimension of the cluster at large distances (up to 
the upper cutoff of the fractal). In addition. F ig. 2b ex- 
hibits a knee- like feature. In the previous work |l3l | a sim- 
ilar knee-like feature served to identify Lifshitz-Slyozov 
length scale L2{t). Here, following Sharon et al. |lj|, we 
subtracted from C{r,t) its initial value C(r, 0), and fol- 
lowed the dynamics of the difference, see Fig. 2c. The 
knee-like feature of Fig. 2b becomes here a local mini- 
mum whose position at different times yield a sharp es- 
timate of L2(t). The "frozen" tail at the distances much 
larger than L2{t), but still much smaller than the system 
size, implies breakdown of DSL 



FIG. 1: Snapshots of coarsening of FVFPs simulated with 
model B. The upper left figure {t = 0) shows a FVFP {D ~ 
1.71) grown in experiment [l^. The rest of the snapshots 
show the simulation results at scaled times t — 290 (upper 
right), 3817 (lower left) and 30000 (lower right). 



Power-law fits of the data shown in Fig 2d yield the fol- 
lowing dynamic exponents: ai ~ 0.24±0.01 for Li ~ t"^ 
and a2 = 0.30±0.01 for L2 ^ t"^ . The same result for ai 
is obtained from a power-law fit of the cluster perimeter 
versus time P{t) - t-°'\ as expected H H [3 III El ■ 
While the value of a2 is very close to that obtained 
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FIG. 2: The dynamics of the equal-time pair correlation func- 
tion C(r, t) (a-c), and the dynamic length scales Li{t) and 
1/2 (t) (d) from simulations with model B. The initial con- 
ditions are the FVFPs grown in experiment [l^. The time 
moments are t = (figures a and b only), 1052, 2950, 13846 
and 3 • 10*. The dashed line describes a power law ~ t^^^ and 
is shown here to guide the eye. See text for further details. 



in earlier simulations of model B |13| and experiment 
with the FVFPs 01 > cn = 0.24 is somewhat larger than 
the values 0.20 — 0.23 reported earlier for these two sys- 
tems IHElIllIillli. Also, noticeable in Fig. 4d is 
curvature of the log-log plot of Li(t). These observations 
put forward a question about the true asymptotic value 
of exponent ai. 

To address this question, we performed a series of 
larger simulations, extending the time interval until t = 
10^ (which is 20 times lon ger than the first phase-field 
simulations of this system 11]). The initial conditions 
for the minority phase were DLA clusters, "reinforced" 
by an addition of peripheral sites. The clusters occupied 
the 1024 X 1024 box; they had a larger fractal range than 
the FVFPs grown in experiment |l4| . Figure 3 shows 
snapshots of the simulated coarsening dynamics. Figure 
4a presents C(r, t) averaged over 6 different realizations 
of DLA. Again, following some initial "evaporation" of 
the minority phase (which happens at an earlier stage of 
the Cahn-Hilliard dynamics), the tail of C{r) is frozen 
until very long times. The dynamic length scale L2{t) is 
shown in Fig. 4b; a power law fit at long times yields 
a — 0.31 — 0.32 which is close to 1/3, as expected. The 
cluster perimeter versus time is shown in Fig. 4c. It can 
be seen that P{t) has not approached yet a power law. 
Therefore, we followed Huse j23| and introduced an ef- 
fective time- dependent exponent —ai{t) which is shown 
in Fig. 4d versus the perimeter P itself. An asymp- 
totic value of —ai{t) is obtained by linear extrapolation 
t ^ 00, that is P ^ 0. This procedure yields ai = 0.34, 
very close to 1/3. 

Overall, our simulations with model B are in remark- 
able agreement with experiment on the Hele-Shaw coars- 
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FIG. 3: Snapshots of coarsening of DLA clusters {D ~ 1.71), 
simulated with model B, at scaled times t — (upper left), 
1350 (upper right), 26591 (lower left) and lO'^ (lower right). 




17345 4 4 4 

10 10 10" 10 10 10^ ^ 1x10 2x10 3x10 



t P 

FIG. 4: C{r,t) at f = 0, 1026, 10521 and lO'^ (a), L2{t) (b) and 
P(t) (c). The dashed line describes a power law ~ t^^^ and 
is shown here to guide the eye. The effective time-dependent 
exponent ai versus P is shown in Fig. d. Linear extrapolation 
to P = (dashed line) yields qi = 0.34. 

ening of FVFPs Breakdown of DSI at large dis- 

tances and the Lifshitz-Slyozov scaling i^/'^ at interme- 
diate distances are firmly established. Our simulations 
show, however, that the "unusual" dynamic exponent 
0.20 — 0.23 is a transient on the way to 1/3. This finding 
explains the apparent independence of the unusual ex- 
ponent on the fractal dimension of the cluster, observed 
in Ref. [l^ . In view of the conjectured universality, we 
expect that the same kind of behavior will be observed 
in a larger-scale experiment on the coarsening of FVFPs, 
and in direct simulations with the unforced Helc-Shaw 
model. 
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